Stability of gap solitons in a Bose-Einstein condensate 
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We analyze the dynamical stability of gap solitons formed in a quasi one-dimensional Bose-Einstein 
condensate in an optical lattice. Using two different numerical methods we show that, under realistic 
assumptions for experimental parameters, a gap soliton is stable only in a truly one-dimensional 
situation. In two and three dimensions resonant transverse excitations lead to dynamical instability. 
The time scale of the decay is numerically calculated and shown to be large compared to the 
characteristic time scale of solitons for realistic physical parameters. 



I. INTRODUCTION 

One of the most fundamental facts about Bose- 
Einstein condensates (BEG) of dilute atomic gases is that 
they can be very well described by the Gross-Pitaevskii 
equation (GPE) 



V'(x,i) (1) 



where -0 denotes the collective wavefunction of con- 
densed atoms with mass M and coupling constant k := 
ATrti^aNA/M , with a being the s-wave scattering length 
and Na denoting the number of atoms in the BEG (see, 
e.g., Ref. 01). Therefore a BEG provides a physical real- 
ization of many nonlinear wave phenomena among which 
the formation of solitons is particularly interesting. 

For our purposes solitons are wavepackets in which the 
dispersive effect of the kinetic term is exactly cancelled by 
the nonlinear interaction energy so that their shape does 
not change. Two fundamental types of solitons have so 
far been experimentally realized in a BEG: dark solitons 
lH, H correspond to a stable density dip in a BEG of 
repulsive atoms. Bright solitons do exist for atoms with 
attractive interaction (k < 0) and arc described by a one- 
dimensional solution of the GPE for vanishing potential 



V'bright(^,<) 



sech(z/w)e 



(2) 



where w = 2h? / {M kid) and uJs — h/{2Mw^). To relate 
the one-dimensional GPE to its three-dimensional ori- 
gin we assumed that the BEG is tightly trapped in the 
transverse direction so that no transverse excitations can 
occur. Taking the BEG to be in the transverse ground 
state amounts to replacing the three-dimensional cou- 
pling constant k by kuj = k/Aj_, where Aj_ ~ 2T:a\ with 
a_L := yjh/ (Mti)j_) is the harmonic oscillator length of the 
transverse potential. The attraction between the atoms 
prevents the dispersion of the sech-shaped wavepacket. 
Note that for kid > 0, corresponding to repulsive atom- 
atom interaction, the kinetic energy and the interaction 
energy have the same sign and therefore cannot cancel 
each other, a bright soliton is then not possible. Very 



recently bright solitons have been created in a quasi one- 
dimensional setup 10, H] where the transverse potential 
V{x, y) tightly confines the BEG, thus suppressing trans- 
verse excitations and three-dimensional collapse. 

In this paper we are concerned with the dynamical 
stability of gap solitons. This collective state, which 
has not yet been realized experimentally, exists for re- 
pulsive atoms {kio > 0) in a periodic potential V{z) 
and is related to bright solitons. The basic idea of a 
gap soliton is the following: as is well known the energy 
eigenvalues of noninteracting particles in a periodic po- 
tential are given by energy bands En{q), where n is the 
band index and q denotes the quasi momentum. If a 
particle's state is prepared in the lowest band only the 
dispersion relation p^/(2M) in free space is replaced by 
the lowest band energy Eo^q). Around the upper band 
edge, which we take to be at g = 0, this energy can 
be approximated by Eo{q) w Eo{0) + q^ /{2M*) where 
M* := [cP Eo{q) / dq^y^\q=o is the effective mass of the 
particle in the periodic potential. Since at the upper 
band edge M* < the "kinetic energy" becomes nega- 
tive and a cancellation with the positive interaction en- 
ergy becomes possible. The corresponding state is called 
a gap soliton. 

Gap solitons have been realized in nonlinear optics us- 
ing a periodic modulation of the propagation medium 
Q. In nonlinear atom optics they have first been pre- 
dicted by Lenz et al. in the context of light-induced 
nonlinearities ^. Here we are concerned with the 
collision-induced nonlinearity appearing in Eq. (^. 

We consider a BEG that is placed in a one-dimensional 
optical lattice, created by far-detuned laser light with 
wavenumber fc^, which produces an optical potential 
of the form V{z) = — Vq cos(2fciz), and is subject to 
a tight harmonic transverse confinement of the form 
Vi^{x,y) — Mlo\{x^ + y^)/2. Here u}±_ is the transverse 
trap frequency and Vo the strength of the optical poten- 
tial. The derivation of the corresponding gap soliton solu- 
tion of the one-dimensional GPE is tedious and includes a 
multiple scales analysis. For nonlinear optics it has been 
derived by Sipe et al. |p^ . For nonlinear atom optics a 
related derivation has been sketched in Refs. |0, |l^, ^ . 
One finds that, within the effective mass approximation. 
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the one-dimensional gap soliton is described by 



III. NUMERICAL METHODS 
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V'gap(2:,t) = ^Phc{z,t)—=sec]\{z/<N)e 



2w 



(3) 



where ip\,c{z^t) = i^bc(^)e~*"^"*-"-'*/'' is the solution of 
the linear Schrodinger equation that corresponds to the 
upper band edge. The parameters w and Cog have 
the same form as w and oj for the bright soliton but 
with M replaced by \M*\ and hid replaced by ki_d := 
I'^iD J \ fhc{z)\'^dz. Apart from (^be the solution (^) just 
corresponds to a bright soliton for a particle of mass \M* \ 
and coupling constant k. 

The range of experimentally promising values for the 
optical potential strength Vb, the number of condensed 
atoms Na, and the transverse confinement frequency uj±_ 
has been examined in a study by Brezger et al. fol- 
lowing standard textbooks [M . It was found that values 
around Na = 400, Vq = h^/{2M) and = 110 s^^ 
should be optimal for the observation of gap solitons. In 
the following we will focus on this case. We will also con- 
sider the effect of a variation of Na which allows to test 
the quasi one-dimensionality of the gap soliton. 



II. STABILITY THEORY OF GAP SOLITONS 
AND BRIGHT SOLITONS 

A stationary solution -00 of the GPE, with chemical 
potential ^, is called dynamically stable if a small devia- 
tion Sij} from ■00 will not grow with time. In this case a 
small perturbation will not cause the solution to evolve 
into a completely different wavepacket. To study dynam- 
ical instability one can either directly integrate the GPE 
or solve the associated Bogoliubov-de Gennes equations 
(EDGE, see, e.g., Ref. |l^). The latter arise when one 
writes the wavefunction V'(x, t) in the form 



^(x, i) = exp(-i/it/?i)(?/;o(x) 6'4>{y.,t)) 



(4) 



and linearizes the GPE in 5^. By making the ansatz of 
a stationary perturbation, 

5V'(x, t) = u(x) exp(— icji) — v*(x) exp(iu;i) , (5) 

one arrives at the BDGE 



Cu — kiPqV 



huju 



-hwv — Cv — k{iPq)'^u 



(6) 



with C := pV(2Af) + V^(x) - ^l + 2k|0oP- A solution 
(m, v) with eigenvalue uj corresponds to a quasi-particle 
mode. The set of all lo forms the quasi-particle spectrum 
which in general is complex. One can show 1 1 1 that if uj 
is in the spectrum then so is lu* . Using Eqn. (5) it is seen 
that the existence of a nonzero imaginary part of one 
quasi-particle frequency implies exponential growth of 
the mode and hence dynamical instability of the state ipQ. 
To demonstrate that the gap soliton is stable therefore 
amounts in showing that the associated quasi-particle 
spectrum is real. 



Although the direct numerical integration of the GPE 
is easy to implement using the split-step method [ p^ it 
has the disadvantage of not being practical for a 3D study 
of the gap soliton. The reason is that the state in- 
cludes two very different spatial scales: the laser wave- 
length 2-K/kL and the width w of the soliton's envelope. 
To cover both scales simultaneously it was necessary to 
consider at least 260 periods of the optical potential or 
2000 spatial points along the z-axis. Since the number 
of total points in 3D is restricted by the capacity of the 
computer a 3D simulation of a gap soliton becomes im- 
practical. We therefore have applied this method only to 
ID and 2D simulations. 

To derive the spectrum of quasi-particles around the 
gap soliton we followed the method of Ref. ||l^ and ex- 
panded the modes (u, v) of Eq. as well as the gap soli- 
ton wavefunction -00 in a- set of basis functions. For the 
z-direction we have chosen a number of Bloch wave- 
functions which are eigenstates of the linear Schrodinger 
equation with potential V{z) = —Vq cos(2fciz). As trans- 
verse modes we used UxUy harmonic oscillator eigenstates 
corresponding to the transverse trapping potential. This 
turns Eq. (|^) into an eigenvalue problem for a 2n x 2n 
matrix, where n :— nxnyUz is the total number of mode 
functions. However, before this eigenvalue problem can 
be solved, one first has to find the exact wavefunction of 
the gap soliton in the given set of basis modes (in practice 
a reasonably large subset is sufficient). 

Since the analytical solution (||) is only approximately 
correct a stability analysis will inevitably lead to a com- 
plex quasi-particle spectrum because of the finite dif- 
ference to the exact solution. To find the exact solu- 
tion '0oxact we have used a self-consistent field approach 
(SCF): after expansion of the GPE in the basis modes 
Eq. (||) is turned into a set of coupled nonlinear algebraic 
equations for the expansion coefficients of ipo- Using the 
approximate solution as an ansatz Atrial. i we insert 
it into the GPE to evaluate the nonlinear terms. The 
resulting equation, 



= (|^ + nx)+«l'/'tnal(x)p)^(x) , 



(7) 



represents a linear eigenvalue problem for -0 and can eas- 
ily be solved using standard numerical methods. We pick 
that solution ?Atriai,2 out of all eigenstates of Eq. which 
has the least deviation from our previous guess V'triai,! 
and iterate this procedure until the change in the trial 
wavefunction is below a given value (we used a relative 
change of 10"^'* as accuracy goal). The converged wave- 
function may then correspond to the true gap soliton. 



IV. ONE-DIMENSIONAL RESULTS 

In one dimension convergence of the SCF algorithm 
can easily be achieved. We have used up to 180 Bloch 
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modes to expand the exact soliton wavefunction which 
covered the upper half of the lowest energy band and the 
lower third of the second Bloch energy band. However 
about 40 modes covering the effective-mass region around 
the upper band edge were sufficient to get about the same 
numerical accuracy. The (real) expansion coefficients for 
the gap soliton can be seen in Fig. |^. It is interesting 
to note that the SCF method essentially amounts to re- 
moving that part of the approximate solution ^ which 
corresponds to the second energy band. The exact nu- 
merical solution therefore is indeed centered around the 
upper band edge of the first energy band where the ef- 
fective mass is approximately constant. We have checked 
whether the final result of the SCF algorithm indeed de- 
scribes a soliton by using it as initial condition for the 
split-step direct integration of the GPE. It was found 
that the solution does not change its shape for times ex- 
ceeding 0.2 s. 
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FIG. 1: The expansion coefficients a, — {ifq(z)\^f}{z)) of the 
wavefunction ^ in the basis of Bloch functions (ji„ around the 
band edge of the first and second band. The circles indicate 
the initial wavefunction and the dots indicate the converged 
wavefunction. The coefficients of the first band are nearly 
identical. The coefficients of the second band are all close to 
zero. The coefficients of the initial wavefunction in the second 
band have been muitiplied by 5 to improve the presentation. 

To verify the stability of the ID gap soliton we have cal- 
culated the quasi-particle spectrum for up to 350 Bloch 
modes as basis functions. The result for the real part 
of the spectrum can be seen in Fig. ^ together with the 
Bloch mode energies. As one can see the two spectra are 
very similar apart from a constant shift. This shift is 
given by the chemical potential and arises because of the 
corresponding phase factor in ansatz (^) The similarity 
of the curves is due to the small number of condensed 
atoms {Na = 400) and the correspondingly small colli- 
sion effects. 

The imaginary part of the quasi-particle spectrum is 
zero except for a single mode (and the corresponding 
mode with complex conjugated frequency) for which the 
frequency is purely imaginary. The value of the imagi- 
nary part for this mode depends on the degree of con- 
vergence of the SCF wavefunction and on the number of 
modes included. For the analytical, non-converged solu- 
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FIG. 2: Real part of the Bogoliubov spectrum and energies 
of the expansion modes for the first two bands in units of the 
recoil energy Erec = h^k\/{2M). 

tion (H) the decay factor is about 0.02 x w^, where wj^ 
is the transverse trap frequency. For a well-converged 
solution it is always very small, below 0.001 x lo^. 

To understand the unstable mode better a few general 
facts about the quasi-particle spectrum are helpful: for 
any potential and any stationary solution of the GPE the 
mode (ii, v) = {ipo, ■0o) is a quasi-particle mode with fre- 
quency oj — 0. This mode is a Goldstone mode associated 
with the symmetry of the energy functional 

^-J{^o(^^ + ^0 + ^^l^ol"} (8) 

with respect to a global phase change ip' = exp(za)V'. If 
the external potential V is absent then, for any stationary 
solution i/jQ of the GPE, there is a second Goldstone mode 
{u, v) — (Vi/'Oi "VV'o)- If arises because of the invariance 
of E against spatial translations. If the effective mass ap- 
proximation would be exact then the gap soliton would 
fulfill the same equation as the bright soliton does in free 
space. Consequently, it would also possess a translational 
Goldstone mode. However, since the periodic optical po- 
tential does explicitly break translational invariance one 
can expect a shift of the complex Goldstone frequency 
away from zero. A second effect that explicitly breaks 
translational invariance is the finite number of Bloch ba- 
sis functions used in the numerical calculations. This is 
equivalent to an optical lattice placed in a box whose 
length is a finite multiple of the lattice period. In our 
case, the box contained 260 periods. 

That the unstable mode indeed corresponds to a mod- 
ified Goldstone mode can also be seen by looking at its 
expansion coefficients (Fig. ||). Obviously the shape of 
u{q) is approximately given by qipo{q) which would de- 
scribe the Goldstone mode if the quasi momentum q is 
replaced by the real momentum as it is done in the effec- 
tive mass approximation. The question remains whether 
this tiny instability results from numerical aberrations, 
from the finite number of periods, or whether it is a real 
physical effect. To shed some light on this question we 
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also have performed a stability analysis of the bright soli- 
ton in free space using the same algorithm. We found a 
similar behaviour: a Goldstone mode develops an imag- 
inary eigenvalue, but this time it is the mode associated 
with the phase transformation. Since one can prove that 
this mode has zero frequency we conclude that the tiny 
imaginary eigenvalue for Goldstone modes is a spurious 
numerical effect and the gap soliton is dynamically stable 
in one dimension. 
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FIG. 3: Expansion coefficients for the modified translational 
Goldstone mode. Shown is the real part (dots) and imaginary 
part (circles) of the function u(z). Numerically it was found 
that v{z) « —u{z). 



V. RESULTS IN TWO AND THREE 
DIMENSIONS 

Having examined the stability of the ID gap soliton it 
is of interest whether it will remain stable in a quasi one- 
dimensional situation. The condition for the latter is usu- 
ally formulated as follows: the interaction energy, which 
leads to a coupling between different modes of the cor- 
responding linear Schrodinger equation, should be much 
smaller than the excitation energy of the transverse trap- 
ping potential. If this is fulfilled a BEG in its ground 
state will effectively behave like a one-dimensional quan- 
tum gas since transverse excitations are off-resonant and 
therefore suppressed. 

This is not the case for a gap soliton, however. The 
reason can be seen by looking at Fig. ^ which displays 
the mode energies of noninteracting atoms in the opti- 
cal lattice and with a tight harmonic transverse confine- 
ment around the upper band edge. The solid line displays 
the energy of atoms in the transverse ground state and 
with longitudinal quasi momentum q in the lowest energy 
band. Each dashed line corresponds to transversally ex- 
cited atoms with a transverse energy of 2huj± to 6hLu±, 
respectively. It is important to observe that there are 
resonances between transversally unexcited atoms with 
q = and transversally excited atoms with q 0. Since 
the gap soliton is a superposition of Bloch modes around 
the upper band edge (g = 0), these resonances have the 



consequence that even for tight transverse confinement a 
true gap soliton does not exist. This situation is qual- 
itatively the same in two and in three dimensions since 
in both cases the free energy levels are given by those 
of Fig. ^. Since in two dimensions the transverse trap- 
ping potential is one-dimensional the multiplicity of the 
energy levels is always one. This is different in three di- 
mensions where a transverse excitation energy of nhujj_ 
has an n -t- 1 fold degeneracy. The number of resonant 
states is therefore larger than in two dimensions. 
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FIG. 4: Energy eigenvalues around the upper band edge for 
noninteracting atoms in an optical lattice and with a trans- 
verse trapping potential. Due to resonances between longi- 
tudinal and transverse excitations the gap soliton will be un- 
stable against transverse decay. The physical parameters are 
given in the text. 

Although a true gap soliton does not exist it is of in- 
terest to examine a quasi gap soliton of the form 

V'quasi = if^gapiz, t)(po{x, y) , (9) 

with ifQ denoting the transverse ground state. It should 
be possible to produce a state like ■(/'quasi using dispersion 
management p^ . Though the quasi gap soliton is not a 
true stationary solution of the GPE, it may be sufficiently 
stable to allow for experimental observation. A signature 
of a quasi gap soliton would be a strongly suppressed 
dispersion of the wavepacket along the z axis. To analyze 
the time scale on which this state decays we first note that 
the transverse ground state has even parity and because 
of parity conservation can only couple to even excited 
levels 2nhjj^. This is the reason why we omitted odd 
transverse excitations in Fig. H. 

We have used the state (|^ as initial condition and 
studied numerically the time evolution of it in two di- 
mensions. To study the influence of the transverse con- 
finement we have considered three BECs with 400, 1600, 
and 25600 atoms in Figs. ^ A, B, and G, respectively. The 
number of atoms Na and the transverse confinement fre- 
quency uj±_ have been simultaneously varied keeping the 
product lo±^Na constant. Gonsequently the interaction 
energy in Eq. (^) is kept constant, since it is propor- 
tional to {ijj^Na)^ ■ This also ensures that the first-order 



soliton condition is fulfilled. The result of the numeri- 
cal time evolution after 23 ms is shown in Fig. ^. While 
some excitations are observable, the state still looks very 
much like a quasi gap soliton. This situation changes af- 
ter 0.23 s (Fig. H). Here one can see strong excitations 
which are growing with decreasing transverse excitation 
frequency. This is a reasonable result since the ratio of 
the interaction energy and the transverse excitation en- 
ergy is greater than one for Fig. ^ B and C hence the 
nonlinear coupling is strong enough to excite the trans- 
verse modes. The figure suggests that in particular the 
transverse modes with 2huj± energy are strongly excited 
because for a fixed value of z there are three density 
maxima in the x-dircction. This is in agreement with the 
qualitative predictions which wc have made above. 




FIG. 5: Density of the wavefunction logjQ \ipf at t = 23ms 
for A) 400 atoms, B) 1600 atoms, and C) 25600 atoms in the 
BEG. See text for details. 

To verify this result we also have calculated the quasi- 
particle spectrum of the state (^) in two and three di- 
mensions. In principle, since the quasi gap soliton is not 
a stationary state, the spectrum is not enough to predict 
the evolution of it accurately and a more sophisticated 
approach is needed [T^ . However, to gain a qualitative 
understanding of the time scale on which V'quasi decays 
the imaginary part of the spectrum is sufficient. 

In 2D we used up to 100 Bloch wavefunctions and up 
to 19 one-dimensional eigenstates of the transverse har- 
monic trap as basis modes. The quasi gap soliton was 
expanded using 51 Bloch states. It turned out that there 
are generally quite many unstable modes, but only few of 
them do have a considerable overlap with the collective 
wavefunction. The number of unstable modes depends 
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FIG. 6: Density of the wavefunction logjQ \tpf at t — 0.23s. 
for A) 400 atoms, B) 1600 atoms, and G) 25600 atoms in the 
BEG. See text for details. 



on the number of basis states used to expand the BDGE 
since the number of resonant transversally excited states 
is growing. However, it turned out that this basis de- 
pendence does only affect modes with a small instabil- 
ity. Some examples are displayed in Fig. 0. Mode ^ 
corresponds very well to the anticipated resonant exci- 
tation of transverse modes. It has a non-zero overlap 
with the quasi gap soliton and otherwise only populates 
even transversely excited basis modes. Correspondingly 
its instability is rather large; only the instability of mode 
, which roughly describes the phase Goldstonc mode, 
decays faster. Mode ^ is unstable because, as described 
above, the quasi gap soliton is not a stationary solution 
of the GPE. However, due to a coupling to transversally 
excited states the decay rate is strongly enhanced as com- 
pared to the non-converged one-dimensional quasi gap 
soliton (the modulus of the transversally excited mode 
coefficients is less than 0.1 and is therefore not visible in 
Fig. 0). Mode is typical for the many unstable modes 
which have no overlap with the quasi gap soliton. There- 
fore, unless transversal excitations are created during the 
experimental preparation of the quasi gap soliton, they 
do not contribute to the decay of it. These modes ap- 
pear because the collective wavefunction provides a lin- 
ear coupling term between transversally excited states 
which are in resonance with each other. Such modes 
can exist for even and odd transverse excitation levels 
without violating parity. Finally, mode ^ describes a 
somewhat off-resonant coupling between the transverse 
ground state and transversally excited states. Because of 
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its off-resonant nature its decay rate Imfo;) is consider- 
ably smaller than for the resonant mode 0A. 



expand the Bogoliubov modes along the z-axis and up to 
30 transverse modes as given in Eq. (|l0|). 
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FIG. 7: Selected unstable Bogoliubov modes for a 2D quasi 
gap soliton. Shown is the energy E (in units of huj±) of the 
basis functions as a function of the quasi momentum q (in 
units of The lowest parabola corresponds to the trans- 
verse ground state. The other parabolas describe a transverse 
excitation of nhu:±,n = 1,2, ... . The thickness of the dots 
corresponds to the modulus of the u or v coefficients with re- 
spect to the corresponding basis function. The smallest dots 
correspond to coefficients of modulus between and 0.1. The 
imaginary part Im(tj) of the unstable modes is given in units 
of the transverse trap frequency uj± = 110 s~^. 
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FIG. 8: Selected unstable Bogoliubov modes for a 3D quasi 
gap soliton. The units are as in Fig. ^, but in addition to 
the basis energy E and quasi momentum q a third quantum 
number m is introduced which is associated with the axial 
angular momentum = mh of the basis functions. 



To analyse the dynamical instability of a 3D quasi 
gap soliton we used again Bloch wavefunctions for the 
longitudinal expansion of the Bogoliubov modes. For 
the two transversal directions we have chosen a basis 
of states which are both eigenstates of the Hamilto- 
nian and the angular momentum operator L^. These 
states can be constructed by using the creation opera- 
tors c± := (aj. ± mp/V2 jl9j. The basis states are then 
given by 



^((n + m)/2)!((n-m)/2)! 



|0> 



(10) 



The energy of these states is given by nhu!± , n — 
0,1,2,... and their angular momentum by mh, m — 
—n, —{n — 2), . . . ,n. We used again 51 Bloch states to 



We found again a large number of unstable modes 
which spuriously depends on the number of transverse 
basis states. However, similarly to the 2D case, modes 
with an instability Im^Lu) > 0.02lli± do not depend sig- 
nificantly on the number of basis states. Some of these 
unstable Bogoliubov modes in 3D are shown in Fig. ||. 
Mode 1^ is again roughly the phase Goldstone mode 
plus some small transversal excitations not visible in the 
figure. Due to the larger number of nearly resonant basis 
states the decay rate is even larger than in the 2D case. 
Mode ||b is again a superposition of resonant states and, 
as in 2D, has the second largest decay rate. These two 
states as well as most of the other unstable modes are 
solely composed of states with zero angular momentum, 
i.e., states which are symmetric in the x-y plane. They 
are therefore very similar to the 2D case. A new kind of 
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instability in 3D occurs in Mode g|C which has no overlap 
with the quasi gap soliton. It is composed out of states 
with m — ±1 and thus describes the decay or increase 
of rotating states linearly coupled by the collective wave- 
function. We found 80 unstable modes of this kind. 



VI. CONCLUSION 

We have examined the dynamical instability of gap 
solitons in a BEC in one, two, and three spatial dimen- 
sions under the condition of tight transversal confine- 
ment. Using different methods we found that a truly 
one-dimensional gap soliton is stable. In higher dimen- 
sions transversal excitations which are resonant with the 
upper band edge forbid the existence of a real gap soli- 
ton. However, a quasi gap soliton may be experimentally 
prepared which behaves like the ID gap soliton for a time 
smaller than the smallest decay time of one of the unsta- 
ble Bogoliubov modes. In 3D we numerically found this 



time to be in the order of 1/(0.133 x uj±) which is the 
decay time of mode A in Fig. ^ For a transverse trap 
frequency of uj± = 110 s^^ we expect the quasi gap soli- 
ton to be sufficiently stable for about 70 ms. This should 
be long enough for experimental observation. 
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